Spin echo without an external permanent magnetic field 
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The spin echo techniques aim at the elimination of the effect of a random magnetic field on the 
spin evolution. These techniques conventionally utlize the application of a permanent field which 
is much stronger than the random one. The strong field, however, may also modify the magnetic 
response of the medium containing the spins, thus altering their "natural" dynamics. We suggest 
an iterative scheme for generating a sequence of pulses which create an echo without an external 
permanent field. The approximation to the ideal echo improves with the sequence length. 

PACS numbers: 76.60.Lz 



I. INTRODUCTION 



The use of echoes is a standard technique in spin reso- 
nance experiments^, where a large number of echo pulse 
sequences are used. Echo sequences typically consist of 
a number of high-frequency pluses that induce controlled 
rotations on the precessing spin. The objective is to re- 
move or reduce dephasing due to inhomogeneities in the 
external magnetic field which give a variation of the pre- 
cession rate across the sample (known as inhomogeneous 
broadening). Since any two-level quantum system can 
be mapped on a spin-^ in a magnetic field, similar oper- 
ations can be performed on any such system, including 
a qubit. For example, there are recent experiments on 
superconducting qubits2*^ (were it is called charge echo, 
since the physical degree of freedom is two charge states 
of a superconducting grain) and electron spins in quan- 
tum dotsi. Common to the existing echo techniques is 
the restriction that in order to apply the proper sequence 
of control rotations one must know the direction of the 
magnetic field. In other words, the sequences are able 
to remove dephasing due to inhomogeneities in the field 
amplitude, but not the field direction (the randomness 
in the direction of the magnetic field is indeed negligible 
if a strong external field is applied). Starting with the 
experiments^ that showed long spin coherence times of 
electrons in bulk GaAs there has been interest in using 
electrons confined in quantum dots as qubits. It appears 
that the limiting decoherence factor in these systems is 
the hyperfine interaction between the electron and nu- 
clear spins&. As long as no special preparation is made, 
the nuclear spin system will give rise to an effective ran- 
dom magnetic field seen by the electron. This field will 
be random both in strength and direction, and in the 
absence of an external strong permanent magnetic field 
the usual echo sequences can not be applied. Because the 
response of the nuclear spin system is much slower than 
the precession of the electron spin we can in a certain 
approximation assume the nuclear field to be constant. 
Motivated by this we will in this paper investigate the 
problem of finding an echo sequence that is applicable 
in the case of a constant effective field with an arbitrary 
and unknown magnitude and direction. The goal of this 



sequence is to return the spin to its starting position at 
the end of the sequence, independently of the magnetic 
field. 



II. FORMULATION OF THE PROBLEM 

Consider a spin-i particle precessing in a constant 
magnetic field. Let n be the unit vector along the pre- 
cession axis and lo be the precession angular frequency. 
The rotation of the spin state during time r is then given 
by the unitary operator 

U — cos \I + i sin x n i a ii (1) 

where I and <Ji are the identity and Pauli matrices, and 
X = lot/2. The 7r rotations about the coordinate axes 
are denoted X — ia x , Y — ia y , and Z — ia z . 

The simplest form of a spin echo is for the situation 
where the spin is precessing in a field with known direc- 
tion, say along the z-axis, but with unknown magnitude. 
Then n x — n y — and n z = 1, while x is arbitrary. The 
usual echo sequence then consists in waiting for the time 
r, applying an X rotation, waiting time r and applying 
a final X rotation. The success of the procedure is ex- 
pressed by the fact that XUXU = —I is the identity (up 
to a global sign, which is unimportant), independently 
of U . By symmetry, the same is true if X is replaced 
by a 7T rotation around any axis in the ccy-plane (per- 
pendicular to the precession axis). In particular we have 
YUYU = -I. 

We want to extend this to the case were n is not known. 
That is, we want to find a sequence of control rotations 
A,B,C,...,F such that 

FU ■ ■ ■ CUBU AU = I (2) 

for any U. We make the following assumptions: i) The 
control pulses are effectively instantaneous, meaning they 
can be performed in a time much less than the preces- 
sion period, ii) The external magnetic field is unknown 
but constant, so that the precession operator U does not 
change in time. Note that we have written Eq. as if 
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the time r between the pulses is fixed. Different time in- 
tervals between the pulses are achieved by choosing some 
of the A, B, C, . . . , F to be the identity. This is sufficient 
if all intervals are integer multiples of a smallest unit. 
Intervals with irrational ratios will require a more gen- 
eral form than Eq. (J2J). In this paper we only use equal 
intervals. 



III. ITERATED MAPPINGS 

As was explained above, if we know the direction n 
of the external field we can create an echo by applying 
two 7r pulses around any axis perpendicular to n. The 
sequence XUXU can be used if n x = 0, and similarly 
we may use YUYU if we know that n y = 0. Consider 
the longer sequence XUYUXUYU for which, using the 
general U of Eq. Q , we get 



XUYUXUYU 
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(3) 



= (8n x n y sin x — 1)1 — 8m x n y sin \ cos x a x 
8in x n y sin 4 xo y - Ain x n y sin 2 x(l - 2n y sin 2 x)er 2 
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We see that XUYUXUYU = -I for either n x = or 
ii y = 0. This may seem like a small gain, but this se- 
quence is the key to the full solution. Let us construct 
the mapping 

U -> T{U) = XUYUXUYU. (4) 

The idea is that T(U), being composed of the arbitrary 
U and the fixed control rotations, will be "less arbitrary" 
than the original U. Iterating this mapping we then con- 
struct the set [/« = T(U), t/ (2) = T(U^), ... of pulse 
sequences that will be better and better approximations 
to the identity. Since contains U^ four times 

the time needed for the sequence U^ is 4™r. That is, 
the length of the sequence grows exponentially in n. This 
is of course unfortunate as the sequences quickly will be- 
come impractically long. However we will see below that 
for a large portion of the space of parameters determin- 
ing U, a few iterations n are sufficient to reach a good 
approximation to the identity. The domain of "bad" pa- 
rameters shrinks exponentially with the increase of n. 

To illustrate this we represent a rotation by the polar 
angle and azimuth angle <f> of the rotation axis (so that 
n x = sin 9 cos <f>, n x — sin 9 sin <fi and n z — cos 9) and the 
rotation angle x- Figure^shows the evolution of the pa- 
rameters (9, <fi, x) upon successive transformations given 
by Eq. (@J. The set of initial rotations fills this space 
uniformly, see Fig. n( a )- Figures H](b)-(f) show the suc- 
cessive iterations U^ - U^ of these points, illustrating 
the convergence of the mapping. 

From Eq. (|3J) we can write down explicit formulas for 
the parameters (n', x') of U' = T(U) in terms of (n, x), 
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FIG. 1: (a) A set of rotations U filling uniformly the space 
(6, (f>, x) °f parameters. Each point corresponds to a certain 
rotation U. (b)-(f) successive iterations of the mapping T(U) 
showing the convergence to the fixed point U = —I. 
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We see that sin% = n x = n v = 0, n z = 1 is a fixed point, 
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and this represents the identity operator. The stability 
of the fixed point can be analyzed by expanding close 
to the fixed point in the independent small quantities 
e s = sin%/\/2, £ x = n x /^/2 and e y = n y /2. We get 
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were the primed quantities refer to the transformed ro- 
tation U'. It is clear that if e s ,e x , e y are small quantities, 
then e' s ,e' x , e' y are even smaller and the fixed point is lo- 
cally stable. 



IV. CONVERGENCE PROPERTIES 

We studied numerically the convergence of the map- 
ping T(U) for points that are not close to the fixed point. 
Looking at Figure ^ we see that although most points 
converge to the vicinity of the fixed point in a few iter- 
ations, there are some points that do not converge fast. 
To get a better understanding we do as follows. Let us 
choose some initial rotation angle x an d for each point 
in the (9,<f>) plane construct the sequence U^, U^ 2 \ ... 
stopping when f/( n ) is within a specified distance from 
the identity (we used the stopping criterion x^ < 10 -3 ). 
Figure El shows the (9, 4>) plane for various x, the shades 
of gray representing the number n of iterations needed 
for convergence. We see that for small x convergence is 
fast and the pattern is simple, but for larger initial x the 
pattern of convergence time is quite complex, and that 
there exist "hard" points, i.e. initial rotations (n, x) that 
need a large number n if iterations to converge. 

To estimate the fraction of initial rotations that need 
a certain number of iterations to converge to the fixed 
point we do the following. We start with an ensemble of 
random initial rotations U characterized by a unit vec- 
tor n distributed uniformly on a sphere and by an angle 
X distributed uniformly from in the inerval [0, it] . We 
proceed with the iterations, Eq. until x^ < 10~ 3 
is reached. The logarithm of the fraction p n of initial 
rotations U that needs n iterations to converge is shown 
in Figure |3 as a function of n. We find that about 74% 
of the initial rotations will converge after applying the 
third iterate and about 91% do so after the iterate 
. Observe that except for the first points, all points 
fall on a straight line which means that the fraction p n 
of "difficult" initial rotations decays exponentially with 
the number of iterations, 



Pn<x.e 



1.1. 



(5) 



If k is the number of repetitions of U in T(U) (for the 
mapping ffl we have k = 4 ), the total time of the se- 
quence TJ( n > is t = rk n , where r is the time between 
the pulses. The relation JSJ can then also be written as 
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(a) X = 0.5. 



(b) X = i-o. 




(c) x = 1-25. 



(d) x = 1-5. 
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FIG. 2: The levels of gray indicate the number n of iterations 

that is needed to reach a rotation with the rotation angle 
x (n) K 1Q -3 The darker shade 

means faster convergence 
according to the scale at the bottom. 



the fraction p(t) of the initial rotations which did not 
converge till time t, 



p(t) oc 
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(6) 



The relations Eqs. iJSJ and JSJ are established by running 
a simulation with the use of a specific convergence cri- 
terion, < 10~ 3 . Changing the criterion affects the 
proportionality coefficients in these relations, but does 
not change the values of exponents a and (3. 

Let M n denote the set of initial rotations that do not 
converge in n iterations. [As an example, the domain Mg 
of "difficult" initial rotations (according to the criterion 
X < 10~ 3 ) is presented in Fig. 0]. The sequence of sets 
M n is such that each set is contained in the previuos one, 
M n C M„_i, and the fraction of points p n in Figurc[3]is 
proportional to the difference in the volumes of M n+ \ and 
M n . The set M = f£L Q M n of "infinitely hard" points 
appears to be a fractal with fractal dimension D« 1.5 
according to the box-counting algorithm we used. 
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FIG. 3: The fraction p n of the initial random set of rotations 
which did not reach the condition x < 10 -3 till the n-th 
iteration. 
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manent field to the system. The advantage of such a 
method is that between the pulses the dynamics of the 
system is not influenced by external perturbations. In the 
context of ESR of a quantum dot, this method may help 
distinguishing between different effects of the hyperfine 
interaction: this interaction creates some random effec- 
tive magnetic field acting on the electron spin, but may 
also lead to electron spin relaxation. The absence of the 
external permanent field in this problem is crucial, as its 
application definitely suppresses the spin relaxation part 
of the hyperfine interaction effects^. 

We have provided a solution to the general echo prob- 
lem in terms of a set of longer and longer pulse sequences 
that give successively better approximations to the iden- 
tity operator. We have tested several longer sequences 
e.g T(U) = ZUXUYUXUZUXUYUXU. In all cases 
we found that the exponent (3 in Eq. @ is independent of 
the particular mapping chosen. Whether this represents 
some inherent property of the problem or only is the case 
for the limited class of mappings we have studied is not 
known. For example, we have not studied sequences with 
control pulses other than tt rotations about the coordi- 
nate axes. We have also not ruled out the possibility of 
a solution that will yield an ideal echo with the help of a 
finite number of control pulses. 

The idea of using iterated mappings to generate pulse 
sequences has been used in NMR applications^ but as far 
as we know this particular problem or the mapping we 
study was never discussed. The mapping we have used 
was also proposed in Ref. in the context of dynamical 
decoupling of a qubit. 



FIG. 4: The set Mg of initial rotations that do not converge 
in nine iterations, according to the criterion \ < 10~ 3 . 
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